All Possible Coupling Schemes in XY Spin Chains for Perfect State Transfer 
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We investigate quantum state transfer in XY spin chains and propose a recursive procedure to 
construct the nonuniform couplings of these chains with arbitrary length to achieve perfect state 
transfer(PST). We show that this method is capable of finding all possible coupling schemes for 
PST. These schemes, without external control fields, only involve preengineered couplings but not 
dynamical control of them, so they can be simply realized experimentally. The analytical solutions 
provide all information for coupling design. 



o 

(N 

a 



Oh' 
I 

=3 : 



> 



PACS numbers: 03.67.Hk, 05.50.+q 



Quantum information and quantum computation can 
process lots of tasks which are intractable with classi- 
cal technologies. Although many schemes such as quan- 
tum dotsflj, ion trapf^l, NMRfsl have been discussed 
extensively, a macroscopic scalable quantum computer 
still seems to need a channel, often known as quantum 
wire, to transmit or exchange quantum states between 
inner parts of the quantum computer. These architec- 
tures require to implement a transmission process for 
an unknown quantum state from one place to another 
which is often called quantum state transfer. In a seminal 
paperQ, Bose proposed a spin chain model, whose evolu- 
tion was governed by a reasonable Hamiltonian, and con- 
sidered the fidelity of state transfer in this model. Similar 
results were also derived by studying dynamical proper- 
ties of^ entanglement transition in Heisenberg XY spin 
chain|5|. This model, in which two processors are con- 
nected through a spin chain as quantum wire, is useful for 
quantum computation based on Heisenberg interactionfil 
or measurements [7]. 

Although some important and significant results have 
been found, see for example [3, all of the available 
results are just concerned with uniform interaction, i.e. 
the couplings between any two nearest-neighbor sites are 
the same. For this case, however, it is shown that when 
N > 4, where N is the number of the sites in XY chain, 
PST is impossible This drawback of uniform inter- 
action motivates people to find some modified models 
to achieve "long" distance PST. Some works considered 
long-range interactions^, and some concentrated on nu- 
merical simulationsllOj . One feasible choice is to preengi- 
neer the couplings[8|, i.e. choose special nonuniform cou- 
plings to achieve PST, and some specific analytical cou- 
pling schemes were found(i,[ill,[i3|. The necessary and 
sufficient conditions for the couplings of PST, which can 
provide a criterion to verify the preengineered schemes as 
well as to find new analytical ones, were derived from a 
more systematical treatment of this problem[13] by mir- 
ror inversion [Til and quantum computation. However, 
all these preengineered schemes are obtained through a 



"verifiable" but not "constructive" way. Thus, we have 
not yet got all possible coupling schemes for PST. 

In this Letter, we start from the necessary and suf- 
ficient conditions of PST. After preselecting the eigeir- 
values of a XY spin chain Hamiltonian, we propose two 
recursive formulas of the couplings for both even and odd 
N cases and prove them by mathematical iirductioir. Fur- 
ther discussions demonstrate that this method is capable 
of finding all possible coupling schemes for PST in XY 
chain with arbitrary length. Experimentally, our PST 
schemes can be realized, for example, by superconduct- 



ing circuits and quantum bus[lj|, nanoelectromechanical 
resonator arraysfl^ or cold-atom optical lattice [l6|. 

Next, we first review some basic concepts of state 
transfer protocol using spin chain as the channel[4i, 1^. 
An unknown qubit, as encoded in site 1, is attached to 
one end of a spin chain when the chain is initialized to 
the all spin-down ground state (state initialization is not 
necessary [l^, and our results can be generalized to these 
cases). Due to the coupling between site 1 and 2, free 
evolution of the system causes the unknown state to dis- 
tribute among the chain. After a specific interval, we 
want to recover this unknown state at the opposite end 
of the chain to achieve state transfer. 

A reasonable Hamiltonian for this task is XY type 
Hamiltonian 
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B^{<7^-1), (1) 



where Ji is the coupling strength between sites i and 
and Bi is the external static potential, or control field, at 
site i. cr^,cr^,cr^ are the three Pauli matrices. One im- 
portant observation is that Hamiltonian ([T]) commutates 
with the total z-spin operator X]i=i '^f ■ Thus, '^f 
a conservation, and the evolution of the system in these 
state transfer cases will just involve the subspace spanned 
by ground state and TV one-site excited states. By the 
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Jordan- Wigner transform, which mapslH]) to 
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XY model can be solved exactly. Hamiltonian(I5]) de- 
scribes an TV-site hopping model subjects to nonuniform 
external fields. Let |i) denotes the single excited state 
at site i, HamiltonianQ in a 2^-dimensional space will 
reduce to an iV-dimensional subspace spanned by 
Explicitly, 
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in {\i)} basis. The fidelity of this transfer procedure can 
be expressed as (iV| e^'-^"^ where r is the time in- 
terval of the free evolution. The equivalent conditions 
for PST, i.e. \{N\e-'"''^ \1)\ = 1, are: (a) the reflec- 
tion symmetry Bi = B^+i-i and Ji = Jn-i- (b) after 

r 

sorting the eigenvalues of Hn— in decreasing order, the 

TT 

difference between any two adjacent eigenvalues is an odd 
number 13|. All schemes discovered before required (a) 
as part of their protocols and designed the eigenvalues of 
Hn- to be {-fc, -fc + 1, • • • fc - 1, fc} for 2k 6 

TT 

{qiP + k) + (2p + l)k} for fc = 0, . . . , or {-k + 

1 
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which all satisfy (b). All these coupling schemes are spe- 
cial solutions for PST, and our main result in this Letter 
is to show how to get all possible couplings for PST in 
the absence of external fields, i.e. Bi — 0. Because of the 
perfect transfer condition (a) and the postulation Bi — 0, 
HaniiltonianQ becomes 
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whose eigenvalues are symmetric about zero. Owing to 
this symmetry, there are only [N/2] independent cou- 
plings and [A^/2] independent eigenvalues inQ (0 is al- 
ways an eigenvalue when N is odd). Our purpose is to 
construct the couplings {Ji} from a set of preselected 
eigenvalues {A^} satisfy (b). We will first consider even 
N cases and show how to derive {Ji} effectively. Then, 
we generalize these results to odd N cases, and finally 
show the completeness of this method, i.e. it can get all 
possible coupling schemes for PST. 



For even TV cases, we assume the eigenvalues of 
are {±Ai, . . . , ±A„} where n = N/2 , A^ G N and Ai > 
A2 > ... > A„ > (if none of {Jj} is zero, then the 
eigenvalues of are nondegenerate 1 1 9j ) . and omit the 

T 

scale factor — . {Ji} and {A^} are connected through the 

TT 

characteristic polynomial of the HamiltonianlU): 



Det(iJjv-A/) = [|(A2-A2) 



(5) 



which, by expanding it with respect to A^, is equivalent 
to a set of equations: 



N-l 



1=1 



ki+i-ki>2 



... T'' — 



2 

3 • ■ 



(6a) 



(6b) 



(6c) 



and we want to derive {J,} from {A;}. This is often 
called an inverse problem. Notice that we still use J^-i 
rather than Ji when i < N/2 despite they are equal 
just for elegance of the expressions. We first introduce 
some notations for convenience. Denote jl^^ — J^/^^ 
and jj^' = (j|^')^ for 1 < i < n — 1 whose mean- 
ing will become clear soon. Here, the superscript N 
denotes the dimension of the matrix and we im- 
ply the eigenvalues of Ho are {±Ai, . . . , ±A£)/2} and 
its couplings are {j]^'}- The main idea is to obtain 
{j'l^^} from when we require the Hamiltonians 

construct by them respectively share the same eigen- 
values {±Ai, . . . , ±A(£)_2)/2}- Further, 

for 1 < k < n. 
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.[N-2] .[N-2] .[N- 
Jn-2 '"Jk-1 



n 



3n-l"'3k + 

- AN -2] 

- J n-l ■ 



where Jq^' = 0, and 
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1 where the products in the 



numerators and denominators involve terms only if the 
indices of them are not larger than n — 1 and n respec- 
tively. With these notations, we will show the following 
equation permits us to get {Ji} from {A^} directly: 
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Eq.lO allows to construct from . • . 

Aj^n \ ■ ■ ■ , jj+j} and A„. Thus, when we know 

by adding one more parameter A„, we can derive 

Jn^'j Jn^'iJ ■ ' ■ jI^' one by one explicitly. Now, we need 
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to prove Eq. ([7]) is consistent with Eq. ^ . Direct calcula- 
tion shows A„ satisfies a continued fraction: 
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1, . . . , n, are all eigenvalues of (U). Fig.([T]) shows idea of 
the proof. 

For odd N cases, we assume A„ > A„_i > . . . > Ai > 
0. Define jlf' = 2{j\^^f and ' = (jf l)^ for 1 < ^ < 
n — l,n— [N — l)/2, and define and as the same 
as even N cases. The corresponding recursive formula 
for odd N is 



Eq.® is equivalent to Det(i?jv — A„/) ~ 0. Actually, by 
expanding Det(i?i — A/) in terms of order i — 1 deter- 
minants, we will find the original continued fraction for 
Det(iJiv - A„/) = is 



+ (-l)^-iA, 



+ (-l)^^A„=0. (9) 



Due to the symmetry between Ji and Jn-ii we can move 
upper half of the continued fraction to the right hand of 
the equal sign. After taking a square root on both sides, 
we obtain Eq.([8]), which means A„ is actually an eigen- 
value of (jlj. The square root operation is exactly the 
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origin of why we denote = j]/^' but 
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for i ^ n before. Next, we will prove Eq.® is correct 
for arbitrary N by mathematical induction. We assume 
the permutations of {Ai, — A2, A3, . . . , (— l)"A„_i} form 
a group keeps { } unchanged which is actually true 

for {jj^'}. The following step is to prove are also 

invariant under the permutation of A„ and —A„_i which, 
with the assumption above, directly induces {±Ai} for 
i — 1, . . . , n are the eigenvalues of Hj^ when { j]^'} are 
constructed from Eq.Q. Obviously, is unchanged 
under the permutation of A„ and — A„_i. For jl^i, 
we can expand it using A„, A„_i and in which 

{ji^ '^^} irrelevant to A„ and A„_i. Notice that this 
expression has similar form with j|j^5^' when j|j^5^' is 
expanded by A„_i, A„_2 and { ^'}, and if we replace 
A„, A„_i and {jj^""*'} in by -A„_i, -A„_2 and 
{■il^ respectively, we will find they are indeed the 
same one. Owing to the assumption that the permu- 
tation of A„_i and — A„_2 keeps j|i^5^' unchanged, we 
conclude that is also unchanged under the permuta- 
tion of A„ and — A„_i. This method, demonstrating the 
invariance by replacement, is applicable for other 
and we can further prove all are invariant under 

the permutation of A„ and — A„_i. 

Combining this proof and the fact that {Ai, — A2} ac- 
tually form a group for {j[^^i2*'}, we can prove the 
permutations of {Ai, — A2, . . . , (— 1)"+^A„} form a group 
keeps unchanged. Furthermore, if A„ is an eigen- 

value of (HI), then, according to Eq.®, {±Ai} for i = 



Am 

Jk 

AN] 



At A 



N 



= A 



N 

fc-1 



fc-1 



(10a) 
(10b) 



where k — n, n — 2, n — 4, • • • till we get The dif- 

ference between n and k in Eq. ()10ap is an even number 
which implies Af^ appears in the right hand of Eq. ([TO|) 
alternately. Just like even N cases, we can directly check 
A„ is an eigenvalue of ^ by the continued fraction rep- 
resentation when {Ji} are expressed by Eq. (|10|) . and the 
factor 2 appears in the definition of jl^^ also comes from 
the continued fraction structure. Although the main idea 
is the same, there are still some differences between even 
and odd cases. First, are no longer unchanged un- 

der the permutation of A„ and — A„_i when N is odd. 
Instead, the permutations of {A„, A„_2, A„_4, . . .} and 
{A„_i, A„_3, A„_5, . . .} form two groups keep in- 
variant respectively (if we consider all the eigenvalues 
{±Ai}, then both even and odd cases have two groups 
formed by interlaced eigenvalues respectively which keep 
{jj^'} unchanged, see Fig.(IT])). Second, it's interesting 



to see 
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structure when a similar replacement as in even N cases 
been made. With the help of this property, we can prove 
A„_i is also an eigenvalues of H^. Combining it with the 
fact that A„ is an eigenvalue of and the symmetry 
property between A„ and A„_2, by means of mathemat- 
ical induction, we assert the eigenvalues of Hn whose 
off-diagonal elements are constructed from Eq. ljlOp are 



- A, - A,- A, -A, ■ • ■ -A |-A OA A , ■•■ A, A, A, A, 
I I ■ , ■■ ■ I I 
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Figure 1: (a) and (b) are recursive procedures for even and 
odd cases respectively. Permutations of {A^} with the same 
colour form one group. Although is always an eigenvalue 
in (b), the permutation group does not contain it. Arrows 
on the straight Unes indicate the recursive directions, e.g. the 
Tightest arrow in (b) impUes if A„ is an eigenvalue of Q when 
{Ji} are constructed from Ea. (|10p then so does A„_2. 
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{0, ±Ai} ioT i = 1, . . . ,n indeed. 

The completeness of this method comes from the fact 
that ^ is uniquely determined by its eigenvalues when 
{Ji} are all positive [2^. This also implies all real cou- 
pling schemes for 1^ are uniquely determined by its 
eigenvalues. Since the completeness is available only if 
all {Ji} are real, we need to prove the positivity of {Jf} 
for Eq.(l7]) and Eq.JTO]). This is more apparent when 
we factor out the common factors of each equation in 
Eq.(l7]) and Eg. pO|) . After factorization, we will find 
each expression contains two factors, one is positive and 
the other is monotone with respect to A„. Considering 
A„_i > A„ > and A„ > A„_i in even and odd cases 
respectively, we assert {jif } are all positive and {Ji} 
are all real which satisfy the completeness condition. In 
a word, Eq.([7|) and Eq.JTO]) are complete for all possible 
coupling schemes. 

Up to now, we have solved both even and odd N cases 
in the absence of external control fields {Bi}. This con- 
structive method allows us to calculate the couplings 
from a set of preselected eigenvalues. We have chosen 
several sets of 50 numbers whose interval between any 
two adjacent ones in each set is a random odd number 
in the domain [1, 100]. In general, we got the couplings 
within 10 seconds. This numerical calculation shows our 
method is effective. Although the resultant couplings of- 
ten have enormous numerators and denominators caused 
by the continued fraction structure of the constructive 
method, we can choose some specific eigenvalues and 
then get compact coupling schemes. For example, choos- 
ing the eigenvalues as {±(T + + ?(2S' + 1))}, where T 
and iS* are two non- negative integers, for i = 1, 2 . . . , ^ 
when N is even, we will find Jf are and 

({l+2T) + a+i)iWS))iil+2T) + (N+l-i)(l+2S)) ^^^^ 

odd i respectively. 

The model used here is also similar to that we en- 
counter in population transfer in an A'^-level system in 
which N discrete energy levels are equivalent to N sin- 
gle excited states {\i)}. Assuming the only interaction to 
be that of electric-dipole transitions and each frequency 
of the laser to be close to resonance with two adjacent 
states, after rotating wave approximation (a general re- 
view of this topic, see[2l|), Hamiltonian of this problem is 
identical to Eq.Q when we treat the dipole interactions 
as the couplings in XY chain. Our results for PST can 
also be used to design the amplitude of each frequency of 
the control laser to achieve perfect population transfer. 

In this Letter, we have considered the problem of trans- 
ferring an unknown state from one end of a spin chain 
to the other end, and proposed two recursive formulas 
for designing the couplings since uniform coupled XY 



chains can not afford PST. We also prove these formu- 
las are complete. Although this method is numerically 
effective, there are still some interesting issues. We set 
the diagonal elements to be zeros, i.e. there is no external 
control field in spin chain or the laser resonances with any 
two adjacent levels in an TV-level system. This is not nec- 
essary for PST or perfect population transfer. Non-zero 
diagonal elements break the symmetry of the spectrum of 
the Hamiltonian, and the eigenvalues no longer appear in 
pairs. Nevertheless, the continued fraction is also avail- 
able when we replace A hy Bi — A. We expect similar 
formula for cases involve control fields which, of course, 
will contain N recursive equations but not [A^/2] for an 
A'^-site spin chain. Another question is whether there are 
other simple coupling schemes for special selected eigen- 
values. We have tested some simple sets of eigenvalues, 
but the couplings still seem complicated. 
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